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Stark shifts of potassium and rubidium Dl lines have been measured with high precision by 
Miller et al In this work, we combine these measurements with our all-order calculations to 
determine the values of the electric-dipole matrix elements for the 4pj — 3d,/ transitions in K and 



•4(^3/2 transitions 



for the 5pj — ^dji transitions in Rb to high precision. The ipi/2 — 3(^/2 and 5pi/2 - 
contribute on the order of 90% to the respective polarizabilities of the np\i% states in K and Rb, and 
the remaining 10% can be accurately calculated using the relativistic all-order method. Therefore, 
the combination of the experimental data and theoretical calculations allows us to determine the 
np — (n — l)d matrix elements and their uncertainties. We compare these values with our all-order 
calculations of the np — (n— l)d matrix elements in K and Rb for a benchmark test of the accuracy 
of the all-order method for transitions involving nd states. Such matrix elements are of special 
interest for many applications, such as determination of "magic" wavelengths in alkali-metal atoms 
for state-insensitive cooling and trapping and determination of blackbody radiation shifts in optical 
frequency standards with ions. 

PACS numbers: 32.70.Cs, 31.15.Ar, 32.10.Dk, 31.15.Dv 



The values of the various electric-dipole matrix ele- 
ments in alkali-metal atoms are needed for the variety of 
applications ranging from reducing decoherence in quan- 
tum logic gates [2| to the study of fundamental symme- 
tries [J IH . The all-order method that includes all sin- 
gle and double (SD) excitations of the Dirac-Fock (DF) 
wave function to all orders of perturbation theory was 
shown to give values for the primary np — ns transitions 
in alkali-metal atoms in excellent agreement with high- 
precision experiments 0. There are many interesting 
applications involving np — n'd transitions in alkali-metal 
atoms and other monovalent systems, but there are very 
few benchmark experiments, such as lifetime measure- 
ments, to compare with theoretical calculations. More- 
over, the only high-precision lifetime measurements of the 
lowest nd states of any alkali-metal atom (g, Q , carried 
out for Cs, are in disagreement [8( with the Stark shift 
values for cesium Dl and D2 lines @, The break- 
down of the correlation correction terms is very different 
for the np — ns and np — n'd transitions, creating a need 
for additional benchmark tests. In this work, we deter- 
mine the values of the <lpj — Zdji transitions in K and 
bp.j — 4dji transitions in Rb to high precision using ex- 
perimental values of the Stark shifts of Dl lines [l| in 
these systems. 

The motivation for this work is two-fold. First, we pro- 
vide the recommended values for these transitions to be 
used for various applications, such as determination of 
the "magic" wavelengths in alkali-metal atoms for state- 
insensitive cooling and trapping and calculation of the nd 
state polarizabilities. Second, we also conduct all-order 
calculations of these transitions in order to carry out a 



benchmark comparison of the accuracy of the all-order 
method. The conclusions reached in this work allow us 
to provide recommended values for a variety of the tran- 
sition properties of monovalent systems and more accu- 
rately evaluate their uncertainties. Such transition prop- 
erties are needed for the evaluation of the black-body 
radiation and quadrupole shifts in ions, light-shifts and 
quadrupole polarizabilities in Ba + which were recently 
measured (ill [l2j, branching ratios for various decay 
channels, and other applications. Such benchmarks are 
also useful for the understanding of the accuracy of the 
all-order calculations conducted for the analysis of the ex- 
perimental studies of parity violation with heavy atoms 
and search for an electron electric-dipole moment. 

The Dl line Stark shifts in alkali- metal atoms were 
measured with high precision by Miller [l[ and Hunter 
et al. @, [IH using a pair of cavity-stabilized diode 
lasers locked to resonance signals. The K and Rb 
measurements, 39.400(5) kHz(kV/cm)- 2 and 61.153 
kHz(kV/cm)~ 2 , respectively, represent a three order of 
magnitude improvement in accuracy upon previous ex- 
perimental results. In this paper, we determine the 
4pi/2 — 3g?3/2 and 5p±/2— 4d 3 / 2 electric-dipole (El) matrix 
elements in K and Rb, respectively, using these experi- 
mental Stark shifts. We also compare these values with 
our all-order calculations of the np — (n — l)d matrix 
elements in K and Rb for a benchmark test of the ac- 
curacy of the all-order method for transitions involving 
nd states. The values of the 4p 3 / 2 — 3dj and 5p3/2 — ^dj 
electric dipole matrix elements in K and Rb, respectively, 
are obtained by combining our recommended values for 
the 4pi/2 — 3^3/2 and 5pi/ 2 — 4e?3/2 transitions with the 
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TABLE I: The contributions to the scalar polarizability for the 4p!/ 2 state in K and their uncertainties. The corresponding 
energy differences [Tj| and the absolute values of the lowest-order (DF) and final all-order electric-dipole reduced matrix 
elements are also listed. The energy differences are given in cm -1 . Electric-dipole matrix elements are given in atomic units 
(eao), and polarizabilities are given in dp, where ao is the Bohr radius. 



Contribution 


k 


{4p 1/2 ||7?||fc) DF 


(4p 1/2 ||73||fc) SD 


E k - E 4pi/2 


«o(4pi/ 2 ) 


a (ns) 


4.s 


4.555 


4.102 


-12985 


-94.8(2) 




5 s 


3.974 


3.885 


8041 


137.3(1.4) 




6s 


0.925 


0.903 


14466 


4.127(3) 




7s 


0.485 


0.477 


17289 


0.962 




8s 


0.319 


0.315 


18780 


0.386 




9s 


0.233 


0.231 


19663 


0.198 




1US 


nisi 


n i so 

U. loU 


zuzzy 


117 














a main (nd 3/2 ) 


3^3/2 


8.596 


7.949 


8552 


540.5(9.7) 




4d 3 / 2 


0.769 


0.097 


14413 


0.05(5) 




5^3/2 


0.105 


0.336 


17201 


0.48(47) 




6^3/2 


0.030 


0.340 


18711 


0.45(30) 




7d 3 / 2 


0.063 


0.296 


19613 


0.33(18) 




8^3/2 


0.069 


0.253 


20193 


0.23(11) 




9d 3 / 2 


0.067 


0.216 


20587 


0.17(8) 


a tail (nd 3/2 ) 










4.5(4.5) 


„ core 

a 










5.5(3) 


Total 










602(11) 



corresponding accurate theoretical ratios. 

We start by expressing the experimental Stark shifts 
as the difference of the ground ns and the first excited 
npi/2 state polarizabilities of the respective atoms [![. 
It is convenient for this purpose to use the system of 
atomic units, a.u., in which e,m c , 4-7ren and the reduced 
Planck constant Ti have the numerical value 1. Polar- 
izability in a.u. has the dimensions of volume, and its 
numerical values presented here are thus measured in 
units of <2q, where ao ~ 0.052918 nm is the Bohr radius. 
The atomic units for a can be converted to SI units via 
a/h [Hz/(V/m) 2 ]=2.48832xl0- 8 a [a.u.], where the con- 
version coefficient is Aire^a^/h and the Planck constant h 
is factored out. 

The Stark shifts in the Dl lines of K and Rb yield 
the following values of the differences of the scalar dipole 
polarizability of the npij 2 and ns states [l[ : 

A K = ao(4pi/ 2 ) - «o(4s) = 317.11(4) 4 (1) 
A Rb = ao(5pi/a) " "o(5s) = 492.20(7) og. (2) 

The static El atomic polarizabilities of the alkali-metal 
atoms are dominated by the valence contribution, a v , 
for all valence states. The small ionic core contribution 
evaluated in Ref. [l6| using random-phase approximation 
approach does not affect the present calculation as it is 
the same for 0:0(^1/2) and ao(ns) and does not con- 
tribute to the Stark shifts. The counter term a vc that 
needs to be introduced to correct ionic core polarizability 
for an occupancy of the valence shell is very small for the 
ns state and is entirely negligible for the np state. The 
valence polarizabilities are given in the sum-over-states 



approach by 



3(2j„ + l)^- E k -E v 



(3) 



where (fc[|D||u) is the reduced electric-dipole matrix ele- 
ment for the transition between states k and v, and Ei 
is the energy corresponding to the level i. The sum over 
the intermediate k states converges very rapidly. In fact, 
the first two ns — npi/2 and ns — np 3 / 2 transitions en- 
tirely dominate the ground state polarizabilities. Since 
these values are known experimentally [T7| . our values 
for the ground state polarizabilities contain very little 
theoretical input. Such a calculation has been described 
before d, [H, [H, H3| and we do not repeat the details 
here. Our resulting values are ao(4s)=289.6(6) Oq and 
an(5s)=318. 35(62) an. These results are in agreement 
with values of [5J, [l8( . The uncertainty comes nearly en- 
tirely from the uncertainties of the experimental values 
of the np lifetimes. 

The polarizabilities of the npi/2 states are dominated 
by a single transition, np\/2 — (n — l)d 3 / 2 , allowing us to 
use Eqs. (JU [2I> to derive the matrix elements that are the 
subject of the present work. As a result, Eqs. ([HE]) can 
be rewritten as 



A K = 



A Rb = 



ao(4Pi/ 2 ) - a (4s) = ^r J - L 



3/2 -~*Pl/2 

o thor ( 4 Pi/2) - 289.6(6) = 317.11(4) (4) 



f. X /pr \ l (4rf 3/2 ||£>||5pi /2 )' 

ao(5pi/ 2 ) - ao(5s) = -— 

3 E 4d 



Ei 



5pi/2 



3 



TABLE II: The contributions to the scalar polarizability for the 5pi/ 2 state in Rb and their uncertainties. The corresponding 
energy differences [J_5| and the absolute values of the lowest-order (DF) and final all-order electric-dipole reduced matrix 
elements are also listed. The energy differences are given in cm -1 . Electric-dipole matrix elements are given in atomic units 
(eao), and polarizabilities are given in ciq, where ag is Bohr radius. 



Contribution 


k 


<5p 1/2 ||D||fc> DF 


(5 Pl /2||D||fc) SD 


E k - E 5pi/2 


ao(5pi/ 2 ) 


main / \ 

a (ns) 


5s 


4.819 


4.231 


-12579 


-I04.ll(l5) 




6s 


4.256 


4.146 


7554 


166.5(2.2) 




7s 


0.981 


0.953 


13733 


A OOf/l/ 1 \ 

4.835(16) 




8s 


0.514 


0.502 


16468 


1.120(7) 




9s 


0.338 


0.331 


17920 


n A A o ( o \ 

0.448(o) 




1 He 


9/17 




lo * oo 


U.ZOU^Z J 




lis 


0.192 


0.189 


19338 


135 (lS 


a tail (ns) 










1.9(0.2) 


a main (nd 3/2 ) 


4d 3 /2 


9.046 


8.017 


6777 


694(30) 




5d 3/2 


0.244 


1.352 


13122 


10.2(9) 




6rf 3 /2 


0.512 


1.067 


16108 


5.2(1.1) 




7rf 3/2 


0.447 


0.787 


17701 


2.6(4) 




8rf 3 /2 


0.369 


0.605 


18643 


1.4(2) 




9d 3/2 


0.307 


0.483 


19243 


0.89(10) 


a tail (nd 3/2 ) 








10.5(10.5) 


Q corc 










9.08(45) 


Total 










805(31) 



+ a# hcr (5p 1/2 ) - 318.35(62) = 492.20(7), (5) 

where we substituted the ground state polarizability val- 
ues and separated the contribution from the np\/2 — 
(n — l)c?3/2 transition. However, the remaining contri- 
butions to the npi/2 polarizabilities grouped together as 
cvg ther (npi/ 2 ) still give 10 % to the polarizabilities of the 
n Pi/2 states and need to be evaluated accurately for our 
approach to yield high-precision values. We describe such 
calculation below. For completeness, we describe the full 
theoretical evaluation of the npi/ 2 polarizabilities first 
and then remove the dominate contribution to determine 

a8 ther (^i/ 2 ) ■ 

We separate the valence polarizabilities into two parts, 
a mam , containing the contributions from the states near 
the valence state, and the remainder a tal1 . We calculate 
the matrix elements contributing to the main term using 
the SD all-order method. We conduct additional semi- 
empirical scaling of our all-order values (SD SC ) where we 
expect scaled values to be more accurate based on the 
analysis of the dominant correlation correction contribu- 
tions. We refer the reader to Refs. [1 S, SH, S3] for the 
description of the all-order method and the scaling proce- 
dure. The experimental energies from [l3|, [l5j are used in 
all main term contributions. The remaining terms from 
highly-excited one-electron states are included in the a tal1 
part. The a tal1 is calculated in DF approximation using a 
complete basis set functions that are linear combination 
of B-splines (23[. We use 70 splines of order 11 for each 
angular momentum state. A spherical cavity radius of 
220 a.u. is chosen to accommodate all valence orbitals 
included in the calculation of a main . We chose to include 
as many states as possible into a mam in order to decrease 



the uncertainty in the remainder term. 

The contributions to the scalar polarizabilities of the 
Apx/2 state in K and 5pi/ 2 state in Rb and their un- 
certainties are listed in Tables Q] and [Til respectively. 
The corresponding experimental energy differences [131 
Il5l . |2J, |25| and the absolute values of the lowest-order 
(DF) and final all-order electric-dipole reduced matrix el- 
ements are also listed. The lowest-order values are listed 
in order to illustrate the size of the correlation correc- 
tion for all transitions. We use the experimental num- 
bers from Ref. [l7| along with their uncertainties for the 
primary ns — npi/2 transitions (for example, 5s — 5pi/2 
transition in Rb) . High-precision relativistic SD or scaled 
SD all-order values are used for all the remaining main 
term transitions. The uncertainties given for the matrix 
elements are equal to the differences between the SD ab 
initio and scaled values. 

The tail contribution is rather small for the ns sum, 
but is significant for the nd^/2 sum. In order to evaluate 
the uncertainties in the tail contributions, we calculated 
a few last main terms by using the DF approximation 
and compared the resulting values with the all-order val- 
ues quoted in Tables UJ and [H] In the case of the ns 
states, the DF values differ from our accurate values by 
only 7-10 % with the difference decreasing with the prin- 
cipal quantum number n. As a result, we assign a 10 
% uncertainty to the ns tail values. In the case of the 
nd^/2 states, the correlation corrections are extremely 
large and nearly cancel the lowest-order contribution for 
K. The differences between the all-order and DF values 
are on the order of 90 % for K and 60 % for Rb. As we 
mentioned above, the large nd 3 / 2 tail contributions are 
the reason for the inclusion of so many states into the 
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main term. We take the uncertainty in the nd^/2 tail 
contributions to be 100 %. 

While such an estimate represents a high bound on 
the value of the Rb tail based on the comparison of the 
DF and all-order values, the case of K requires some 
additional consideration owing to larger discrepancies 
of the DF and all-order values even for n = 9. We 
have conducted additional all-order calculations of the 
K ~ El matrix elements with n > 9. We 

found that our calculational scheme becomes impracti- 
cal for n > 15. Such a problem is expected because the 
lowest-order energies of the higher basis set states are 
large and positive resulting in cancellations in the de- 
nominators of the MBPT terms and consequently lead- 
ing to the divergence of the all-order iteration procedure. 
In our approach, the tail does not sufficiently converge 
at n — 15 to significantly reduce the tail error. As a re- 
sult, we chose a different approach to ensure that we do 
not significantly underestimate the tail in the K calcula- 
tion. We compared our final results for K and Rb with 
experimental values for the ropi/2 state polarizabilitics 
by combining the Dl line Stark shifts from [l| with rec- 
ommended ground state polarizability values from [l8j |. 
Our K and Rb values differ with the experimental re- 
sults 606.7(6) a.u. and 810.6(6) a.u. by 0.69 % and 
0.77 %, respectively. If the tail in the K calculation were 
significantly underestimated, we should have seen a sig- 
nificant mismatch of the K and Rb comparison since the 
tail problem is specific to K calculation. 

Tables U and [Til show that the uncertainties of the 
ndz/2 tail values give the overwhelmingly dominant con- 
tributions to the uncertainties of cvq jher (np 1/2) • 

Subtracting the contribution from the np 1 / 2 — (n — 
1)^3/2 states from our final theoretical values for npi/ 2 
polarizabilitics and removing the corresponding uncer- 
tainties from the total error budget, we obtain 



, other 



, other 



(4p 1/2 ) = 61.6(4.8) al 
(5p 1/2 ) = 111(11)4 



Substituting the ag thcr values and experimental energies 
from [24( into Eq. ((4]), we obtain the following absolute 
values of the El matrix elements: 



K (4p 1/2 ||£>jj3d 3/2 ) 
Rb (5p 1/2 ||X>||4d3 /2 > 



7.984(35) ea Q 
8.051(63) ea . 



(6) 
(7) 



The uncertainties of these values are essentially defined 
by the uncertainties of the respective nd^/2 tail contribu- 
tions for both K and Rb. The contributions of the uncer- 
tainties in the Stark shift measurements and ground state 
polarizabilities to the uncertainty of the recommended 
matrix elements values given by ([6l [7]) are negligible. 

We compare these final recommended results with our 
ab initio and scaled values in Table IIIII Since the con- 
tributions of the triple excitation are important for the 
accurate evaluation of these matrix elements, we also con- 
ducted another all-order calculation partially including 



TABLE III: Comparison of the recommended values for the 
n Pi/2 — (n — 1)^3/2 electric-dipole reduced matrix element 
in Rb and K derived from Stark shifts in this work (listed in 
row labeled "Present work" ) , with our theoretical calculations 
carried out using single-double all-order method (SD), single- 
double all-order method including partial triple excitations 
(SDpT), and their scaled values. Absolute values are given. 
Units: eao. 

K(4p 1/2 - 3d 3/2 ) Rb(5p 1/2 - 4d 3/2 ) 

Present work 7.984(35) 8.051(63) 
SD 7.868 7.846 

SD SC 7.949 8.017 

SDpT 7.956 7.994 

SDpT sc 7.953 8.015 

Final theory 7.949(80) 8.02(17) 



TABLE IV: Recommended absolute values of the np — (n— l)d 
electric-dipole reduced matrix element in K and Rb. Units: 
eao. 



K 



Rb 



4pi/ 2 

4 P3/2 

4p: 



3/2 



3d 
3d. 
3d 



3/2 
3/2 
5/2 



7.984(35) 
3.580(16) 
10.741(47) 



5??l/2 
5J53/2 



4d, 



4d 



3/2 
3/2 



8.051(63) 
3.633(28) 



5p 3 /2-4d 5/2 10.899(86) 



the triple excitations to the extent described in Ref. [f|. 
We refer to these results as SDpT values in Table |TTTJ and 
text below. The corresponding scaled values are listed 
in row labeled "SDpT sc " . We take the SDpT sc v alues as 
out final values (see, for example, Refs. |22j, [26( for the 
discussion of this issue). We note that SD SC and ab initio 
SDpT values are essentially the same in the case of K and 
very close in the case of Rb. The uncertainty of the final 
values is taken to be the maximum difference between the 
final values and SD, SDpT, and SDpT sc ones. Our all- 
order values are in excellent agreement with the values 
derived from the Dl line Stark shift. We also conclude 
that our procedure for the uncertainty evaluation of the 
theoretical matrix element values for the np — (n — l)d 
transitions overestimates the uncertainty, especially for 
Rb. 

We also evaluated the recommended values of the 



4P3/2 — 3dj El matrix elements in K and 5p 3 / 2 — 
Adj El matrix elements in Rb using our recom- 
mended values ([6l [7]) and appropriate theoretical ra- 
tios. The ratios (4 Pl ^\\D\\2,d i/2 ) / {Ap m \\D\\Zd m ) and 



<4p3 /2 ||£>||3d 3/2 / 



The ratios (4pi/ 2 | 

2 } / (4^3 /2 1 1 D 1 1 3^5/2) inK are essentially in- 
dependent of the correlation correction, i.e. the DF 
and all-order values are the same to better than 0.1%. 
Therefore, the theoretical values of the ratio are accu- 
rate enough for such a recalculation. The case of Rb is 
exactly the same as that of K. The complete set of our 
recommended values for all six np — (n — l)d transitions 
considered in this work is given in Table HVl 

In summary, relativistic all-order calculations of the 
nsi/2 and npi/2 state polarizability are presented. The 
calculations are combined with the experimental Stark 
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shift values in order to determine the np±/2 — (n— l)c?3/2 
matrix elements in K and Rb atoms with high precision. 
The values of the matrix elements calculated using the 
experimental Stark shifts are found to be in excellent 
agreement with the values of the reduced matrix elements 
evaluated using the relativistic all-order method. This 



work provides a benchmark test for the all-order matrix 
elements involving nd states of monovalent systems. 

This research was performed under the sponsorship of 
the National Institute of Standards and Technology, US 
Department of Commerce. 
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